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Abstract 

A nonlinear description for perpendicular particle diffusion in strong electromag- 
netic fluctuations is developed by using the fundamental Newton-Lorentz equation. 
Although not based on the same approach, the recently presented nonlinear guid- 
ing center (NLGC) theory is recovered as a special case. The approach used here is 
rather based on the argument that a well defined particle gyromotion does not exist 
in strong fluctuations than on the assumption that the particle gyrocenter follows 
magnetic field lines which themselves separate diffusively. The assumption of a guid- 
ing center motion and the diffusive separation of magnetic field lines is absolutely 
central to the NLGC theory. It is argued that the NLGC result should provide most 
accurate results for strong fluctuations. Furthermore, as a direct consequence of the 
particle equation of motion, it is shown that particle diffusion in one perpendic- 
ular direction is governed by the fluctuations in the other normal direction. This 
results contradicts the NLGC result, where perpendicular diffusion is triggered by 
fluctuations in the same direction. This is of particular interest for anisotropic per- 
pendicular diffusion in a non-axisymmetric turbulence. Future numerical simulation 
results for non-axisymmetric magnetic turbulence and their comparison with the ap- 
proach presented here and the NLGC theory have to provide an answer whether 
particle diffusion in one perpendicular direction is governed by the fluctuations in 
the other normal direction or by fluctuations in the same direction. 
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1 Introduction and Motivation 



Diffusive particle motion in random fields plays a key role in space physics, 
astrophysics and the physics of fusion devices. Turbulence properties are es- 
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sential for understanding the three-dimensional (anisotropic) diffusive particle 
transport in a collisionless, turbulent and magnetized plasma such as the so- 
lar wind or the interstellar medium. Of particular interest are the transport 
coefficients describing particle diffusion perpendicular to an ambient magnetic 
field. 



In spite of its long-standing importance in space and astrophysics, perpendic- 
ular diffusion has been an unsolved puzzle during many decades and a variety 
of studies have been carried out to achieve closure and to pin it down at a 
theoretical level. Mo dels hay e been proposed for hard-sphere scattering in a 
magnetized plasma ( Gleeson . 19691 ) and for relat ed extensions developed on 
the basis of the Boltzmann equation ( Jonesl . 1990f) . but they seem to be inap- 
plicable to spa ce plasmas where (electro)mag netic fluctuations trigger particle 



scattering (see Bieber and Matthaeusl . 1997 ) 



Other models are based on the field line random walk (FLRW ) limit emerg- 
ing, for slab turbul e nce g eometry from quasilinear theory fQLT Uokipii . 1966; 
Jokipii and Parkerl 19691) . which has been considered and used in subsequen t 
studies fe.g. lForman et al.l . ll974lFormarJ . ll977tlBieber and Matthaeusl . Il997). 
Although FLRW provides for a physically appealing picture, it has been shown 
recently that its applicability is questionable for particle transport in certain 
turbulence geometries, particularly th ose with at least on e ignorable coordi- 



nate being the case for slab geometry ( Jokipii et al. . 19931 ). Furthermore 



nu- 



merical simulations, taking into account the magnetic nature of the turbulence, 
have show n that FLRW fails to explain perpendicular tran sport of low-energy 



nave snown tnat f L.KW tails to explain perpendicular transport ot low-energy 
particles ( Giacalone and Jokipii . 1999L iMace et al. . 2000 ) , since pitch-angle 



scattering is neglected in the FLRW limit (e.g.'lQ in et al. 



2002al ). 



It has been argued for some time that diffusion along a background magnetic 
field can reduce perpendicular diffusion to subdiffusive levels if the t urbu- 
lence reveals slab geometry only (see, e.g., IUrchl . ll977trK6ta and Jokipii l200tl 



Qin et all Eo02bl ). However, when the turbulent magnetic field has sufficient 



structure normal to the mean magnetic field, subdiffusion as seen in pure sla b 



turbulence can be overcome and diffusion is recovered (see iQin et all l2002af) 



Recently, iMatthaeus et all J2003) proposed the so-called nonlinear guiding 
center (NLGC) theory. For their approach, they use the following two key 
assumptions: First, perpendicular transport is governed by the velocity of 
particle gyrocenters that follow magnetic field lines. Second, the magnetic 
field lines themselves separate diffusively. Based on these two assumptions, 
Matthaeus et al. obtain for the velocity of the particle gyrocenter in the re- 
direction the expression 



5B* 
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where 5B X is the x-component of the turbulent magnetic field, Bq is the 
strength of the background magnetic field. The constant a has to be deter- 
mined after the fact. Consequently, in the NLGC theory, particle diffusion in 
normal direction is solely triggered by the diffusive separation of the underly- 
ing magnetic field line s. Employing the so-cal l ed Taylor-Green-Kubo ( TGK 



formulation (see, e.g., Bieber and Matthaeusl . Il997t iMatthaeus et all EoO 



for more details), Matthaeus et al. derive, on the basis of equation (1), the fol- 
lowing nonlinear integral equation for the perpendicular diffusion coefficient: 



a 2 v 2 

w 2 



d 3 k 



v/Xz 



c k\ + K zz k 2 z + 7 (k) ' 



(2) 



Here, k zz = v\ zz /3 is the parallel diffusion coefficient, v is the particle speed 
and X zz is the parallel mean free path. The wavevector components parallel 
and perpendicular to the mean magnetic field are given by k z and k\ = k^ + k 2 , 
respectively, and S xx is the spectral amplitude of magnetic fluctuations in the 
x-direction. The decorrelation rate 7(k) allows to include dynami cal effects 
due to the decay of turbulent energy (see, e.g.. iBieber et al. . 1994 . for more 
details concerning dynamical magnetic turbulence). Since the NLGC theory is 
based on the TGK formalism, equation (2) only admits Markovian diffusion; 
anomalous transport processes such as subdiffusi on and superdiffus ion can not 
be accommodated by the NLGC model 
discussion). 



see also lZank et al.L I2004L for further 



Consi dering the particular case of a static (7(k) = 0) turbulence. IMatthaeus et al 
f)2003[) use the NLGC approach for weak (5B 2 /B 2 = 0.04) and stronger 
(5B 2 /B 2 = 1.0) turbulence and compare the results with numerical simu- 
lations. The numerical simulation results used for the comparison are ob- 
tained by applying a fourth-order Runge-Kutta method with sub sequent in 



tegrat io n to the equation of mot ion of a charged test particle (see iQin et al. 



Matthaeus et al 



2003). The finding of the comparison is that the 
NLGC theory seems to be consistent with the numerical simulations for the 
case a = l/v^3- Since the NLGC approach has been fruitful for a better un- 
derstanding of perpendicular diffusi on, it was recently applied to cosmic ray 



modulation and shock acceleration (|Zank et all 12004). 



Although the NLGC approach seems to be in agreement with the numeri- 
cal simulation results, several questions arise. Among those are the following 
two: (1) Does a well defined particle gyromotion and, therefore, a gyrocenter 
actually exists for fluctuations being that strong that the background mag- 
netic field has no control over the particle motion? Wandering in the particle 
pitch-angle and gyrophase are then very rapid, and the concept of a gyrocenter 
velocity following magnetic field lines is expected to fail for strong fluctuations. 
(2) By using the fundamental particle equation of motion itself, is it actually 
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possible to derive spatial diffusion coefficients being in agreement with the 
NLGC result? A closer inspection of the ansatz (1) shows that the NLGC ap- 
proach is not developed on the basis of the Newton-Lorentz equation, which 
describes particle diffusion in turbulent electromagnetic fields. The NLGC re- 
sult, equation (2), is rather based on the diffusivity of the underlying magnetic 
field lines. The diffusivity of m agnetic field lines is given by the magnetic fiel d 
streamline equation (see, e.g., Ijokipii and Parkerl . Il969l iRuffolo et ail 120041) . 



which apparently was employed to result in equation (1). 

In view of the two questions given above, a nonlinear description for per- 
pendicular diffusion in strong electromagnetic fluctuations is presented which 
explicitly excludes an "ordered" particle gyromotion. The approach used here 
is based on the particle equation of motion and contradicts the above men- 
tioned key assumptions made by Matthaeus et al. for the NLGC theory. 



2 Basic Equations and Formulation 



For a random, diffusive particle motion and statistically homogeneous and 
stationary flu ctuations, k xx is given by the Taylor-Green-Kubo (TGK) formula 
( Kubol . : 957) which can be expressed as 



x / ^mok(O), 



(3) 



where v x is the x-component of the particle velocity. An analogous expression is 
given for the diffusion coefficient in the ^-direction, K yy . For a large turbulence 
coherence time £, the second-order velocity correlation function (v x (0)v*(£)) 
must go to zero, and the integral in equation (3) approaches a constant value 
for t — > oo. This is characteristic for Markovian diffusion. 

The use of the TGK formula is somewhat critical for situations where particle 
transport reveals rather an anomalous diffusion process than normal Marko- 
vian diffusion. In this case, the standard Fokker-Planck coefficient (anomalous 
transport law) k xx = (Ax 2 ) /2t applies, where the mean square displacement 
scales more general as (Arc 2 ) oc t s . Depending on the exponent s, different 
diffusion processes can be taken into account: s — 1 for the diffusive regime, 
corresponding to a Gaussian random walk of particles; s = 2 in the superdif- 
fusive reg i me, i .e. strictly scatter-free propagation of the particles (see, e.g. , 
McCombl. Il990h. and s < 1 (e.g. 8=1/2, as shown bv bin et all (Ho02b) 



and Kota and Jok oiil §000)) in the case of particle trapping (subdiffusion or 



compound diffusion). For s = 1 (normal Markovian diffusion), the standard 
Fokker-Planck coefficient is equivalent to the TGK formula (3) for the large t 
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limit. 



2. 1 Equations of Motion 



Generally, the random walk of a particle with mass m, charge q and velocity 
v(t) = x(t) in fluctuating fields is governed by its equation of motion, i.e. the 
Newton-Lorentz equation 

p = i 7mv (0] = gE[x(t),i] + ^v(t) x B[x(f),t], (4) 
at c 

where 7 is the Lorentz factor. The magnetic field B = B + 5B consists of 
a uniform, steady background magnetic field B = B e z to which is added 
a broad band spectrum of magnetic fluctuations 5B[x(i),£]. Analogously, the 
electric field E = Eo + 5E consists of a vanishing background electric field 
< E >= E = with superimposed electric fluctuations <5E[x(t),t]. Equation 
(4) can formally be integrated in time, to obtain for the perpendicular velocity 
components the expressions 



t 



«x(*) = — / dt' (c5E x [x(t'),t'\ +v y (t')6B z [x(t'),t'} -v z (t')5B y [x(t'),t'\) 
7771c J 

t 

+ n j dt'Vyit') (5) 





t 

Q 



= ^ / dt'{cSE y [x(t'),t'} +v z {t')6B x [x(t'),t'} -v x (t')6B z [x(t'),t'}) 







-n j dt'v x (t'), (6) 


where Q = qB^j^mc is the relativistic gyrofrequency. 

The first integrals in equations (5) and (6) describe the random velocity com- 
ponents of the particle in the x and 77-direction, respectively, due to the fluc- 
tuating fields. They enter equations (5) and (6) solely by these contributions 
and represent a stochastic force acting on the particle and forcing it to de- 
viate from its unperturbed orbit given by the second integrals. In a uniform 
and ordered background magnetic field, each unperturbed orbit is a perfect 
helix along which the particle moves at constant speed v z , while the transverse 
velocity components v x and v y oscillate at angular frequency Q around a mag- 
netic line of force. Now consider the effects of fluctuations on the unperturbed 
particle orbit. For weak (5B 2 /Bq <ti 1) turbulence conditions, the fluctuations 
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cause the particle to wander in its pitch-angle and gyrophase. The force in- 
duced by B dominates the particle transport and stochastic contributions 
resulting from the first integrals are comparably small. Consequently, a parti- 
cle "remembers" its helical orbit and performs a diffusive but ordered motion 
along a preferred direction defined by Bo- The deviation from its helical orbit 
in perpendicular directions is then marginal, and diffusion is dominated by 
parallel scattering. 

In the opposite case of strong (5B 2 /B 2 > 1) turbulence, the ordered force 
induced by B is small and particle transport is dominated by the "chaotic" 
contributions given by the first terms in equations (5) and (6). Then, a particle 
moves barely on a perfect trajectory and "forgets" its ideal orbit very quickly. 
Decorrelation from the helical orbit then occurs very rapidly, the wandering in 
pitch-angle and gyrophase is then purely chaotic, and a well defined particle 
gyromotion does not exist. Whether the first or second term in equations (5) 
and (6) dominate depends then solely on the magnitude of the ratio 5B 2 /B 2 . 

The decisive step taken here is to assume strong turbulence conditions, that 
the particle has no memory of its formerly helical trajectory. To do so, the 
contribution representing the helical orbit is dropped. This does not imply 
that Bo = 0, it is rather assumed that the random velocity components (first 
integral) overwhelm the unperturbed particle orbit (second integral). Pro- 
ceeding, the irregular fields are expressed by Fourier transforms and the first 
integrals are, for simplicity and the reason of overview, approximated by their 
integrands. This yields 

v x (t) = [ rf^k-xw [ c ^( k) t) + v ( t )5B z (k, t) - v z (t)SB y (k, £)] , (7) 

7mc J 

v y {t) = [ d 3 ke 1 ^ [cSE y (k,t)+v z (t)5B x (k,t)-v x (t)SB z (k,t)] .(8) 
7mc J 



Here, t x and r y denote timescales related to the deviation of the particle po- 
sition from the unperturbed orbit in the x and ^/-direction, respectively. The 
timescales t x and r y might depend on particle rigidity and turbulence proper- 
ties, but they are here introduced for the reason of dimensionality. A future 
study requires, of course, in any case a detailed treatment of the time integra- 
tions. This is far beyond of the scope of the current paper, since the evaluation 
is a study on its own. In particular, the TGK formula (3) requires the eval- 
uation of v x at time zero, which is a result of the time translation invariance 
property usually employed for the deriv ation of the TGK formula (see, e.g., 
McCombl . 1199(1 iKota and Jokkril HoOO) . Obviously, the substitution of equa- 



tions (5) and (6) into the TGK formula would then yield zero for t = 0. 
Actually, the approach used here by employing directly the particle equation 
of motion requires the more general Fokker-Planck definition k xx = (Ax 2 )/2t, 
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and this is currently under investigation. The purpose of this paper is rather 
to demonstrate that the same mathematical structure on diffusion coefficients 
as those of the NLGC model can be obtained if the Newton-Lorentz equation 
is used for the assumption of electromagnetic fluctuations overwhelming ef- 
fects related to the background magnetic field. As becomes clear below, the 
timescales r x and r y are then related to the "numerical" constant a given in 
equation (1). 



2.2 Velocity Correlation Functions 



The substitution of equations (7) and (8) into the TGK formula (3) converts 
the integrands into a sum of nine individual contributions involving random 
particle velocity components and/or turbulent field components at two differ- 
ent wavevectors. Since the particle velocity components are random quantities, 
i.e. (v i) = 0, four contributions vanish. Proceeding, it is assumed that the par- 
ticle velo city components are unc orrected with the local magnetic field vector 
(see also Matthaeus et al. . 2003f ). Furthermore, supposing that the probabil- 
ity distribution of the particle displacements and that of the Eulerian velocity 
field are statistically i ndependent, which is also referred to as Corrsin ' s inde - 
pendence hypothesis (Corrsin, 19591 : Matthaeus et al. . 20031 McCombI Il990). 
one obtains for the velocity correlation function in x-direction the expression 



Mo)<(0> 



-^-V [d 3 k [ d 3 k' (e-*' *® 
■fine J J J x 



c R X x(k, k , £) 



+(v z (o)v:(0)s yy (k,k',o + ( Vy (o)v;(0)s zz (Kk',o 
-(v y (o)v* z (0)s zy (k,k',o - (v z (o)K(0)s yz (k,k',o 



(9) 



and an analogous expression for (v y (0)v Here, the spectral amplitudes 



S nm (k,k r ,0 = (5B n (k,0)5B* m (k',O) 



(10) 



R nm (k,k',C) = (SE n (k,0)5E* m (k',0) 



were introduced, where the subscripts refer to the Cartesian coordinates. Fur- 
thermore, without loss of generality, the initial condition x(0) = was used. 
For statistically homogeneous turbulence, it is convenient to employ the rela- 
tions 5 nm (k, k', f) = 5(k-k')S nm (k,£) and R nm (k, k', ^) = 5(k-k')R nm (k, £), 
and it is assumed that electromagnetic fluctuations can be represented by a 
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superposition of N individual plasma wave modes, i.e. 



A? 



3=1 



(12) 



JV 



5E(k,t) = ^5E J '(k)exp(-^t). 



(13) 



Here, oUjQs.) = ^-^(k) + ^-(k) is a complex dispersion relation of wave mode 
j, where Uj t R is the real frequency of the mode. The imaginary part, r 3 -(k) < 0, 
represents dissipation of turbulent energy due to plasma wave damping. Note 
that such an an s atz for the time variation is usually used in quasilinear theory 
see, e.g., iGarvl |l993; ISchlickeiserl . 12002] ) . where fluctuations are supposed to 
be very weak. Here, however, it is assumed that the same approach holds 
for fluctuations dominating the influence of the background magnetic field on 
particle motion. 



Restricting the considerations to transverse fluctuations, i.e. 5W • k = 0, and 
using Faraday's law, the turbulent electric field can easily be expressed by the 
corresponding magnetic counterparts, yielding 



(14) 



and, subsequently, the relation 



c 2 k* 



(15) 



For simplicity, it is assumed that S 3 zy (k) = S ] yz {k) = for following calcula- 
tions. In what follows, one arrives at 



M0K(O> 



qr x 



\2 jy 

) J2 / d 3 ^e~ jk ' x( ^e l ^ H(k)e+r ' (k)? 
/ i=i 



x 



CO j 



k A 



+M0K(O>^(k) + ( Vy (ox(0)si z (k) 



(16) 



and obtains an analogous expression for (v y (0)v *(£))• 
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2.3 Nonlinear Diffusion Coefficients 



, x as well as K yy . To do so, 



One can now proceed to derive expressions for k 
the time-integration in equation (3) has to be performed, and further progress 
requires the variations of the individual contributions in the coherence time 
£. First, the velocity correlation functions in equation (16) can be written as 



see also Matthaeus et al. . 2003) 



MO)<(0> = (v 2 /3)exp(-</Kn) 



(17) 



where K nn = v\ nn /3 is used, being in agreement with the TGK definition (3). 
Here, X nn is the mean free path in the direction of one of the Cartesian coordi- 
nates. Second, the ensemble-averaged exponential expression can be expressed 
as 



oo oo oo 

(exp [— zk • x(£)]) = J dx J dy J dzP(x,y, z)e 



where the distribution function P(x, y, z) represents displacements of the par- 
ticle position from the unperturbed orbit in Cartesian space due to the decorre- 
lation of the trajectory. The decisive step is to assume a Gaussian distribution, 
i.e. 



P(x,y,z) 



(27r) 3 / 2 a x a y a. 



exp 



x 

2a 2 



2a 2 



(19) 



where the variances a 2 



{Ax 2 ), a 2 y = (Ay 2 ) and a 2 z 
to be diffusive in the sense that (Ax 2 ) = 2k xx £, (Ay 2 ) 
2k zz £. Upon substituting the Gaussian distribution function (19) into equation 
(18), one arrives, after elementary integrations, at 



(Az 2 ) are assumed 
2k w ^ and (Az 2 ) = 



(exp [-zk ■ x(f )]) = exjp[-(n xx kl + K yy k 2 + K zz k 2 z )£\. 



(20) 



Using these two entries, the ^-integration in equation (3) is elementary. Upon 
extending in equation (16) the factor in front of the wavevector integral by 
-Bo (not assumed to be zero), one finally obtains for k 
expressions: 



, x and Kyy the following 



2 2 N 
a x V 



35o 2 U 



d 3 k 



3M 2 
v 2 k* 



A 



A2 + <R 



(21) 
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a 2 y v 2 A 
K yy ~ or2 

^0 j=l 



(A + v/\ zz y + ul R zzK ' (A + v/\ yy f + uj\ r _ 

iV . r o| . 12 A 

^I d3k i^ (W k ) + #&( k )) i^t^ (22) 



Here, A = k^A^ + + « ZZ A^ ~ Tj(k) and a x = r^fi as well as a y = r y Q are 
introduced. Equations (21) and (22), the main results of this paper, form a set 
of two nonlinear and coupled equations for the two diffusion coefficients in the 
two distinguished perpendicular directions. Each diffusion coefficient depends 
not only on k zz , but also on the perpendicular transport parameter describing 
(diffusive) particle transport in the other normal direction. 

An interpretation of the individual contributions can be given by the "right- 
hand" rule known from electrodynamics. The scattering of a particle in a 
certain Cartesian direction can be considered as the result of a stochastic 
force acting on the particle and forcing it to move diffusively in this direction. 
For instance, consider the contributions appearing in k xx . According to the 
Lorentz- force term in equation (4), a "chaotic" force in x-direction requires a 
random particle speed along the z-axis. In the second contribution in equation 
(21), we see that there is a stochastic velocity component v z present, through 
X zz . This random speed, together with the irregular magnetic fluctuation in the 
7/-direction, via S yy , results in the diffusion in the x-direction. Particle diffusion 
in one perpendicular direction is, therefore, governed by the fluctuations in the 
other normal direction. The third contribution results from the turbulent field 
component aligned along the weak background magnetic field. According to 
the Lorentz-force term, a random particle speed in ^/-direction is then required 
or, equivalently, the mean free path X yy . Finally, the first term in equation (21) 
results from Faraday's law by which the fluctuating electric field in x-direction 
is expressed by the corresponding magnetic field components. 

It is instructive to compare the gener al diffusion coefficient (21) with the 
NLGC result by Matthaeus et al. ( 2003f ) in equation (2). For their derivation, 



they assume that the fluctuations are purely magnetic. The limit of vanishing 
random electric fields can be achieved by choosing a^^k) = in equation (21), 
initially derived for the plasma wave viewpoint. The first term in equation (21) 
results from Faraday's law and, therefore, would not occur for purely magnetic 
fluctuations. Consequently, it is dropped for the comparison. Since the concept 
of a superposition of individual wave modes does not apply, N = 1, and 
the corresponding j = 1-sub (super) scripts are omitted. For the dynamical 
behavior of purely magnetic fluctuations, Matthaeus et al. introduced the 
quantity 7(k). To take this into account for the comparison, the plasma wave 
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dissipation rate Tj(k) < has to be replaced by — 7(k). For purely magnetic 
fluctuations, one then obtains the general relations 



' </</ 



a 2 x v 2 



SB 2 



d 3 k 



d 3 k 



Syy(k) 



+ 



S zz (k) 



_B + v/X zz + 1 (k) B + v/X yy + 7 (k) 
B + v/X zz + 7 (k) B + v/A^ + 7 (k) 



(23) 
(24) 



where B = n xx k 2 + n yy k 2 + K zz k 2 . For their approach, Matthaeus et al. neglect 
fluctuations along the background magnetic field, implying here S zz = 0. Their 
assumption of an axisymmetric turbulence can be achieved by using t x = r y , 
S X x = Syy and k xx = K yy . The diffusion coefficient (23) then reduces to the 
NLGC result. The only difference is the term in front of the integral, and 
a comparison leads to the important result that a = t x Q. Matthaeus et al. 
argue that their factor a is related to the gyrocenter velocity. Numerically, 
they found a = l/v^3, implying here t x VL = 1/V3. But this is characteristic of 
strong turbulence, where an "ordered" and well defined gyromotion of charged 
particles is not expected to be present anymore. Equations (21) and (22) and 
their simplified versions (23) and (24), respectively, are, therefore, referred 
to as the strong nonlinear (SNL) theory. The considerations presented above 
and the formal agreement of the simplified SNL equations with the NLGC 
result indicate that the latter will probably provide for most accurate results 
if strong fluctuations are assumed. 



3 Summary and Conclusions 



The focus of this paper is nonlinear perpendicular diffusion of charged parti- 
cles in irregular electromagnetic fields. Using the fundamental Newton-Lorentz 
equation for the random motion of a test particle and assuming that a well 
defined particle gyromotion does not exist in strong turbulence, nonlinear in- 
tegral equations are derived representing particle diffusion in the two distinct 
perpendicular directions. The two new nonlinear diffusion coefficients are valid 
for a plasma wave dispersion relation, i.e. real frequency of a plasma wave 
mode and its decay due to dissipation, depending arbitrarily on wavevector. 
Both transport parameters are not only coupled with the diffusion coefficient 
for parallel particle transport, but they are also coupled among one another. 
The new set of nonlinear perpendicular diffusion coefficient s is referred to as 



the st rong nonlinear (SNL) theory, and the NLGC theory bv lMatthaeus et al. 



(2003) is recovered for the case of purely magnetic fluctuations. Although 
the SNL theory generalizes the NLGC approach and, furthermore, is based 
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on another physically appealing approach, namely the Newton-Lorentz equa- 
tion, a complete theory requires a prescription for determining the input pa- 
rameters a x = t x Q and a y = r y Q or, better, the detailed evaluation of the 
time-integrations in equations (5) and (6). Of particular interest is the obser- 
vation that particle diffusion in one perpendicular direction is governed by the 
fluctuations in the other normal direction, e.g., 

^XX 

S yy . This is a direct con- 
sequence of the Lorentz-term appearing in the Newton-Lorentz equation (4). 
On the other hand, the NLGC result, equation (2), predicts k xx oc S xx . This 
difference might be of interest for perpendicular diffusion in non-axisymmetri c 



turbulence, since S xx and S yy are then not equal (e.g.. iRuffolo et all 120011) . 
Future numerical simulation results for non-axisymmetric magnetic turbulence 
and their comparison with equations (23) and (24) as well as the NLGC theory 
should then provide for an answer whether particle diffusion in one perpen- 
dicular direction is governed by the fluctuations in the other normal direction 
or by fluctuations in the same direction. 
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